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At high magnetic field, the semiclassical approximation 
which underlies the Ginzburg-Landau theory of the mixed 
state of type II superconductors breaks down. In a quasi-lD 
superconductor with an open Fermi surface, a high magnetic 
field stabilizes a cascade of superconducting phases which 
ends in a strong reentrance of the superconducting phase. 
From a microscopic mean-field model, we determine the ther- 
modynamics and the excitation spectrum of these quantum 
superconducting phases. 

PACS numbers: 74.20-z, 74.70-Kn, 74.90+n, 74.60-w 



It has recently been proposed that quasi-one- 
dimensional superconductors should exhibit an unusual 
phase diagram in a high magnetic field [QJD . In a quasi- 
lD conductor (weakly coupled chains system) with an 
open Fermi surface, the magnetic field does not quan- 
tize the semiclassical orbits which are open but it in- 
duces a dimensional crossover in the sense that it tends 
to confine the wave functions along the chains. This 
quantum effect of the field strongly modifies the phase 
diagram predicted by the Ginzburg-Landau-Abrikosov- 
Gor'kov (GLAG) theory which is based on the semiclas- 
sical phase integral (or eikonal) approximation Q. The 
usual Ginzburg-Landau (GL) regime is followed, when 
the field is increased, by a cascade of superconducting 
phases separated by first order transitions which ends in 
a strong reentrance of the superconducting phase where 
the chains interact by Josephson coupling. The supercon- 
ducting state evolves from an Abrikosov vortex lattice in 
the GL regime towards a Josephson vortex lattice in the 
reentrant phase. Between these two limits, the ampli- 
tude of the order parameter and the current distribution 
show a symmetry of a laminar type while the vortices still 
describe a triangular lattice. The cascade of phase transi- 
tions originates in commensurability effects between the 
periodicity of the order parameter and the crystalline lat- 
tice spacing. This high-field-superconductivity can sur- 
vive even in the presence of Pauli pair breaking because 
the quasi- ID Fermi surface allows one to construct a 
Larkin-Ovchinnikov-Fulde-Ferrell state which can exist 
far above the Pauli limited field. A very important as- 
pect is that the temperature and magnetic field scales 
are determined by the coupling between chains [Q. This 
means that the temperature and magnetic field ranges 
where high-field-superconductivity is expected can be ex- 
perimentally accessible if appropriate (i.e. sufficiently 



anisotropic) materials are chosen. 

In this letter, we derive the thermodynamics and the 
excitation spectrum of these quantum superconducting 
phases from a microscopic model in the mean-field ap- 
proximation. We consider a strongly anisotropic super- 
conductor described by the dispersion law (h = ks = 1 
throughout the paper) E(k) = v(\k x \ — kp) + t z cos(k z c) 
where the Fermi energy is chosen as the origin of the en- 
ergies, v is the Fermi velocity for the motion along the 
chains and t z is the coupling between chains separated 
by the distance c. For a linearized dispersion law, the 
y direction parallel to the magnetic field does not play 
any role (as long as the Cooper pairs are formed with 
states of opposite momenta in this direction) so that we 
restrict ourselves to a 2D model. We assume that the 
zero field critical tempeature T c q is smaller than t z so 
that the superconductor is really 3D: in the GLAG de- 
scription, there is no Josephson coupling between chains 
even at T = 0. We consider singlet pairing but we neglect 
the Zeeman term (i.e. we put the g factor equal to zero) . 
The Pauli pair breaking effect can easily be incorporated 
in the present description. 

In the gauge A(Hz, 0, 0), the one-particle Hamiltonian 
obtained from the Peierls substitution Ho = E(k — > 
— iV — eA) is given by 

v(—iad x — kp) + amuj c + t z cos(—icd z ) , (1) 



n 



where a = + (— ) labels the right (left) sheet of the Fermi 
surface and fa is the (discrete) position operator in the 
z direction. We have introduced the energy lu c = Gv 
where G = —eHc is a magnetic wave vector and H the 
external magnetic field perpendicular to the system. 7Yo 
can easily be diagonalized by noting that the momentum 
along the chains is a good quantum number and by taking 
the Fourier transform with respect to m. The eigenstates 
and the corresponding eigenenergies are 

Cj(r) = ^=e*^_ m (at), (2) 



-k x ,l,i 




(3) 



where r = {x,m), t = t z /uj c , a = sgn(fe x ) and J; is the 
Ztli order Bessel function. L x is the length of the system 
in the x direction. The spectrum consists in a discrete set 
of ID spectra. The state , is localized around the Zth 
chain with a spatial extension in the z direction of the or- 
der of tc which corresponds to the amplitude of the semi- 
classical orbits H . Note that the states <^ ; can be ob- 
tained from the localized states introduced by Yakovenko 
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Pj by a gauge transformation. The superconducting in- 
stability can be qualitatively understood from the spec- 
trum (|^). In zero-field, time-reversal symmetry ensures 
that -E'i(k) = E±(— k) so that the pairing at zero total 
momentum present the usual (Cooper) logarithmic singu- 
larity which results in an instability of the metallic state 
at a finite temperature T c q. A finite magnetic field breaks 
down time-reversal symmetry. Nonetheless, we still have 

e fex h T = e 9~J-fc* h I f° r 1' x = — (h+h)G. Thus, whatever 
the value of the field, some pairing channels will present 
the Cooper singularity ~ ln(27f2/7rT) (7 ~ 1.781 and 
is the cutoff energy of the attractive interaction) if the to- 
tal momentum along the chain q x is a multiple of G. This 
results in logarithmic divergences at low temperature in 
the linearized gap equation which destabilize the metallic 
state at a temperature < T c < T c0 Besides the 

most singular channels which present the Cooper singu- 
larity, there exist less singular channels with singularities 
~ In |f2/nu; c | (n ^ 0) for T <C ui c . In this high field limit 
(lu c 3> T), a natural approximation consists in retaining 
only the most singular channels. Such an approximation 
has been used previously in the mean-field theory of the 
field-induced-spin-density-wave phases observed in some 
organic conductors (Single Gap Approximation) || and 
in the mean-field theory of isotropic superconductors in 
a high magnetic field where it is known as the Quantum 
Limit Approximation (QLA) [||. In the following, we 
shall adopt this latter designation. 

In order to obtain the critical temperature when uj c 3> 
T, we consider the two-particle vertex function in the rep- 
resentation of the eigenstates of T~Lq. In the chosen gauge, 
the total momentum q x along the chains is a constant of 
motion. Since in the QLA we consider only the most sin- 
gular pairing channels, the center of gravity of the Cooper 
pair in the perpendicular direction is related to the to- 
tal momentum q x by L12 = (l\ + Z2) /2 = —q x /(2G) as 
explained above, and therefore becomes also a constant 
of motion. Thus the two-particle vertex function can be 
written as x, (^12^12) with q x (L) = —2LG, where 

I12 = l\ — li and 1' 12 = l[ — 1*2 describe the relative motion 
of the pair in the z direction. In the ladder approxima- 
tion, the integral equation for T then reduces to (/ = li 2 , 

I' = I'n) 

+ \x{0) E ^r'^#u(lV), (4) 

a",l" 

where Vf" = a 1 a' 1 Vi t v and 
/" 27r dx 

Vi,i>= — J l (2tcosx)J l ,(2tcosx). (5) 
Jo 2n 

— XV"" is the matrix element of the local electron- 
electron interaction — XS(ri — r 2 ) (A > 0) in the rep- 



resentation of the states z . Note that VJ"" is in- 
dependent of the center of gravity of the Cooper pair. 
x(0) = N(0) ln(27f2/7rT) is the pair susceptibility at zero 
total momentum in zero field and N(0) is the density of 
states per spin at the Fermi level. The preceding integral 
equation is solved by introducing the orthogonal trans- 
formation Uiji which diagonalizes the matrix Vij> . One 
obtains: 

where V is the diagonal matrix U^VXJ . The metal- 
lic state becomes instable when a pole appears in the 
two-particle vertex function which leads to the critical 
temperature 

Tc = ^ e AN<0) %,'o , (7) 

where Vi j is the highest eigenvalue of the matrix V. 
The critical temperatures are shown in Fig^ for the two 
highest eigenvalues of V. This clearly indicates that there 
are two lines of instability competing with each other and 
leading to a cascade of first order transitions in agreement 
with the exact mean-field calculation of T c Q]. Except 
for the last phase, T c calculated in the QLA is several 
orders of magnitude below the exact critical tempera- 
ture: it has been pointed out previously that the QLA 
strongly underestimates the critical temperature |7j] . The 
existence of two lines of instability results from the fact 
that Vij> = if I and I' do not have the same parity. 
Diagonalizing the matrix Vim is then equivalent to sep- 
arately diagonalizing the matrices Pai,2!' and V»+:L,2i'+i- 
In the following, we label these two lines by Iq — 0, 1 so 
that Vi j Q = max; V 2 i + i ,2i+i - Since 2L = h + l 2 and 
h — h have the same parity, L is integer (half-integer) for 
Iq = (Iq = 1) and can be written as L = —Iq/2 +p with 
p integer. Correspondingly, we have q x {L) = (Iq — 2p)G. 
It is clear that the instability line Iq corresponds to the 
instability line Q = IqG which was previously obtained in 
another approach 0] . From Eq. (^) , one can see that the 
superconducting condensation in the channel q x (L), L, Iq 
corresponds to the following spatial dependence for the 
order parameter 

a,l 

(8) 

where we note a = —a. Noting that the matrices V 
and U have a range of the order of t (i.e. Vu', t/y are 
important for |Z|, |Z'| < t), one can see that A & (£),£,io has 
the form of a strip extended in the direction of the chains 
and localized in the perpendicular direction on a length 
of the order of ct. This is not surprising since ^g x iL),L,lo 
results from pairing between the localized states t . 
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Following the original approach proposed by Abrikosov 
||, we construct the order parameter for T < T c 
as a linear combination of the solutions (||): A(r) = 
J2l l(^)^q x (L),L,i { r ) where 2L must have the parity of 



In. Since A 



an extension of the order of ct, a natural choice for the 
coefficients -f(L) is to take 'j(L) ^ if L = —ln/2+pN' (p 
integer) where the unknown integer N' is expected to be 
of the order of t. In order to correctly describe the trian- 
gular Josephson vortex lattice in the last phase (t <C 1) 
pL we choose j(L) = j(p) — 1 (i) for p even (odd) which 
leads to (noting N = 2N') 

i, P 

xJ p N.+ lo _ m (-t), (9) 

where the amplitude A is chosen real. Eq.(^) defines a 
variational order parameter where the two unknown pa- 
rameters A and N have to be determined by minimizing 
the free energy. It can be seen that |A(r)| has period- 
icity a x = 2n/NG and a z = Nc so that the unit cell 
contains two flux quanta: Ha x a z = 20 o (when a tri- 
angular lattice is described with a square unit cell, the 
unit cell contains two flux quanta) . In Ref. , the order 
parameter was constructed by imposing that it describe 
both the triangular Abrikosov vortex lattice in weak field 
(lu c <C T) and the triangular Josephson vortex lattice in 
very strong field (oj c ^S> t z ). Both approaches lead to the 
same order parameter when only the Cooper singularities 
are retained. 

In order to derive the thermodynamics and the excita- 
tion spectrum in the superconducting phases, it is nec- 
essary to determine the normal and anomalous Green's 
functions from the Gor'kov equations 

(iu - W£)G«(r, r', u>) - A CT (r)^(r, r', u>) = S(r - r') , 
Ha, - Hf)Fp (r, r» + A* (r)G? (r, r', w) = , (10) 

with the 
self-consistency equation A*(r) = AT^^i^^r, r, u). 

Here A CT (r) = A (r)V#(r)) where the C( r )'s are 

fcrmionic operators for particles moving on the sheet a 
of the Fermi surface. Aj(r) = — A;(r) is the variational 
order parameter defined by (||). G" and are the 
Fourier transforms with respect to the imaginary time r 
of the correlation functions — (T T tf)" (r, r)ip" ' (r,0)) and 

-(T r ^" t ( r , T )'C t (r,0))- In ©, it is assumed that the 
magnetization M = (B— H)/47r is equal to zero. This ap- 
proximation is justified in the last phase (t <C 1) where M 
is of the order of t 2 and its contribution to the Gibbs free 
energy G(T, H) of the order of t A . In the other phases, we 
expect that the approximation B = H will give reliable 
results at least not too far from the reentrant phase. 

In order to have a simple description of the supercon- 
ducting state, we introduce the magnetic Bloch states 



-Y> 



'q m +pNG,l-pN 



(11) 



where N z is the number of chains. The eigeneneries 



is localized in the z direction with are e 



q is restricted to the first mag- 
netic Brillouin zone: \q x \ — kp G] — Tr/a xl Tr/a x ] and 
q z G] — ir/a z ,n/a z ]. There are N branches {—N/2 < 
I < N/2) at the Fermi level. The order parameter (^) is 
entirely described by the pairing between the states <^ l 
and 0^ G -q -i -i where G = (G, 0). Consequently, the 
Gor'kov equations ( flo|) become diagonal in the represen- 
tation of the magnetic Bloch states and their solutions 
are 



G?(q,M 



e^ ff 2 + |A ? (q,0| 2 ' 
A?(q,0* 



C// + |A«(q,0l 



2 ' 



(12) 



where the pairing amplitude A"(q, I) is defined by 
A«(q,Z) = J ^r^W^^^^W^A^r) 

= A a a l °Vi 0t i *^2"f P e~ ipq * a *U2i+i +pN,i ■ (13) 



Here A 



-A 



The functions G" and 



F^ appearing in (12) are the Fourier transforms 
of the correlation functions — (T T b^ t CT (T)6q \ CT (0)) and 



-(Trb^rK^-^M where b^ a (b^) are 
annihilation (creation) operators of a particle with spin 
a in the state l . 

From the knowledge of the Green's functions, we can 
calculate the free energy of the system. Close to T c , the 
difference between the free energies of the superconduct- 
ing and normal states can be obtained in a GL expansion 
F N [A] = AA 2 + BA A /2 where 



A = 



B 



S 

-Y 



f d 2 v 


A(r) 




A 



T ^ |A«(q,0/A| 5 



3? 



|Ag(q,Q/A| 4 

2N2 



(14) 



c q,Z,er 



where S = L x N z c is the area of the system. The two pre- 
ceding equations can be further simplified by using (fl3|). 
Minimizing the free energy i*jv[A] with respect to A, we 
obtain Fjy = —A 2 /2B. In the reentrant phase where the 
approximation B = H is justified, we find that the min- 
imum of Fpf is obtained for N = 2 || . When the field 
is decreased from its value in the reentrant phase, the 
system undergoes a first order phase transition and the 
minimum of F^ is then obtained for N = 4. This re- 
sult is in agreement with Ref. O where it is argued that 
the first order phase transitions are due to commensura- 
bility effects between the crystalline lattice spacing and 
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the periodicity of the order parameter. Unlike what was 
expected (2|, the best value of TV switches to 6 before 
reaching the next first order transition. This indicates 
the importance of the screening of the external field in 
calculating the free energy in the phases TV > 4. 

The specific heat jump at the transition is obtained 
from AC = —Td 2 F/dT 2 . The ratio AC/C N where C N 
is the specific heat of the normal state is found to be 
always smaller than the (BCS) zero-field value and dis- 
continuous at the first order phase transitions (FigJ^). 
The discontinuity can be related to the slope AT/ AH of 
the first order transition line || : the slope is positive for 
the last transition and negative for the other transitions. 
The magnetization is obtained from M = —dF^/dH. 
Since T c is determined by Fn = 0, M has the sign of 
dT c /dH Jl(J. Each phase will therefore first be param- 
agnetic and then diamagnetic for increasing field, except 
the reentrant phase which is always paramagnetic JTlJ] . 

From the Green's functions (|l^), we deduce the 
quasi-particle excitation spectrum E™ l a = ±(cq j + 
|A£(q, 0| 2 ) 1/2 - A gap 2A£(q, I) opens at the Fermi level 
in each branch I. The spectrum is shown in Fig.|] for 
the last three phases TV = 2, TV = 4 and TV = 6. In 
the very high field limit (t <C 1), the spectrum is almost 
flat, the dispersion of the quasi-particle band being of the 
order of P. When the field is decreased within a given 
phase, the dispersion increases. The minimum excitation 
energy decreases when TV increases. Thus FigJ| clearly 
shows how the system evolves from a quasi-lD (quasi-2D 
if the magnetic field direction is taken into account) be- 
havior in very high magnetic field (i « 1) towards the 
GL regime (u) c <C T) where the spectrum is known to be 
gapless [|. 

In conclusion, we have solved the BCS theory for a 
quasi-lD superconductor in a high magnetic field. The 
theory can easily be extended to include the pairing chan- 
nels which are not considered in the QLA: the results 
presented in this letter are not qualitatively modified. 
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FIG. 1. Solid lines: critical temperature vs magnetic field 
for lo = and lo = 1 in the QLA. Dashed line: exact 
mean-field critical temperature. 



FIG. 2. Ratio r = (AC/C n )/(AC/C n )bcs vs magnetic 
field. 

FIG. 3. Excitation spectrum in the phases TV = 2, TV — 4 
and TV — 6. The units are chosen so that max \A"(q z , Z)| = 1. 
In the phase TV, there are TV/2 distinct branches. 



[1] A.G. Lebed', JETP Lett. 44, 114 (1986); L.I. Burlachkov, 
L.P. Gor'kov and A.G. Lebed', EuroPhys. Lett. 4, 941 
(1987). 

[2] N. Dupuis, G. Montambaux and C.A.R. Sa de Melo, 
Phys. Rev. Lett. 70, 2613 (1993); N. Dupuis and G. Mon- 
tambaux, to appear in Phys. Rev. B. 

[3] See, for example, Superconductivity, edited by R.D. Parks 
(Dekker, New- York, 1969). 

[4] V.M. Yakovenko, Sov. Phys. JETP 66, 355 (1987). 

[5] K. Maki, Phys. Rev. B 33, 4826 (1986). 



4 



